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Abstract In this paper we analyse the relations between a previously described oblate Jaffe model for an ellip- 
soidal galaxy and the observed quantities for NGC 2974, and obtain the length and velocity scales for a relevant 
elliptical galaxy model. We then derive the finite total mass of the model from these scales, and finally find a 
good fit of an isotropic oblate Jaffe model by using the Gauss-Hermite fit parameters and the observed ellipticity 
of the galaxy NGC 2974. The model is also used to predict the total luminous mass of NGC 2974, assuming 
that the influence of dark matter in this galaxy on the image, ellipticity and Gauss-Hermite fit parameters of this 
galaxy is negligible within the central region, of radius O.SRq. 
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1 Introduction 

Galaxies are usually modelled in terms of their morphology, surface luminosity profile, mass-light ratio, velocity 
dispersion profile and the line-of-sight velocity distribution curve (often termed the rotation curve). Recently, 
dynamical modelling of the ellipsoidal galaxy NGC 2974 was performed by Cinzano and van der Marel (1994), 
based on the Gauss-Hermite fit parameters, i.e. the rotation velocities v, velocity dispersions a and deviations 
of the velocity profiles from a Gaussian, as quantified by the Gauss-Hermite moments hj and /14. A similar 
modelling of the ellipsoidal galaxy M32 was given by Qian et al. (1995). After that, Baes and Dejonghe (2004) 
provided a completely analytical class of dynamical models for spherical galaxies and bulges with central black 
holes. Baes et al. (2005) also studied the dynamical structure of isotropic spherical galaxies with a central black 
hole. Some dynamical models of luminous and dark matter for 17 early-type galaxies were recently presented by 
Thomas et al. (2007). Previously a number of more extensive studies have been published (e.g. Gerhard 1993a,b; 
van der Marel & Franx 1993; Bender et al. 1994; Gerhard et al. 1998; Kronawitter et al. 1999; Saglia et al. 2000; 
Halliday et al. 2001). In particular, the works of Emsellem, Goudfrooij and Ferruit (2003) and Krajnovic et 
al. (2005) are very thorough investigations of NGC 2974. The former shows a two-arm gaseous spiral in the 
inner 200 pc of the early-type galaxy NGC 2974 by fitting the TIGER and long-slit stellar kinematic data of the 
galaxy with a dynamical model, not requiring the addition of a disc or a central dark mass. For their dynamical 
model, a numerical two-integral distribution function was derived from a multi-Gaussian expansion mass model 
by using the Hunter-Qian (1993) algorithm. The paper of Krajnovic et al. presents a very thorough analysis of 
NGC 2974 using completely different data and methodology. Krajnovic et al. first obtained a multi-Gaussian 
expansion mass density, fitting simultaneously the ground-based /-band image and dust-corrected PC part of 
the WFPC2/F814 photometry of this galaxy; they then constructed a three-integral model from the recovery of 
the mass density using the Schwarzschild (1979) orbit-superposition method; finally, they provided a dynamical 
modelling of stars and gas in this galaxy by use of the Gauss-Hermite fit parameters. In this paper we again 
model dynamically NGC 2974, by employing distribution functions from self-consistent density-potential pairs 
of an oblate Jaffe model (Jiang 2000) and using the Gauss-Hermite fit parameters and the ellipticity (Sandage 
and Tammann 1981). The line-of-sight velocity distribution curves are calculated from the distribution functions 
for the oblate Jaffe model and quantified by its Gauss-Hermite fit parameters. The derived Gauss-Hermite fit 



Z. Jiang 

Department of Matiiematics, Zhongshan University, Guangzhou 510275, China 
E-mail: mcsjzl@mail.sysu.edu.cn 



2 



parameters for the oblate Jaffe model correspond very well to those of NGC 2974. The model is also applied 
to predict the total luminous mass of NGC 2974, since in the region of a radius of O.SRe {Re = 33.7 arcsec) the 
dynamical influence of the dark halo of NGC 2974 is negligible (cf. Cinzano and van der Marel 1994). 

As indicated above, a distribution function can be derived from a mass density by either the Hunter-Qian 
algorithm or the Schwarzschild method. As in some other published studies (e.g. Qian et al. 1995; Emsellem, 
Goudfrooij and Ferruit 2003), we obtain below a distribution function using the Hunter-Qian (1993) algorithm. 
It is worth mentioning that the first step of our dynamical model is to give directly the oblate Jaffe mass model by 
using only the galactic ellipticity in order to mimic the image of NGC 2974. This is different from the approach 
used by all the earlier works mentioned above. 

NGC 2974 has been studied frequently (e.g. Gallouet et al. 1975; Sandage and Tammann 1981; Demoulin- 
Ulrich et al. 1984; Forman et al. 1985; Canizares et al. 1987; Davies et al. 1987; de Vauouleurs et al. 1991). 
It appears to be completely isolated (Kim et al. 1988), is classified as E4 by Sandage and Tammann (1981), 
and has a regular distribution of gas (Kim et al. 1988) and a dark halo (Cinzano and van der Marel 1994). The 
total mass of NGC 2974 within a galacto-centric radius of 13 kpc is estimated as 3.8 x 10^' Ii^^Mq by Kim 
et al. (1988), and its total luminous mass as 1.2 x lO'^ H^^Mq (Cinzano and van der Marel 1994), where h is 
Hubble's constant. A range of a factor of 2 is given for h [h is in units of 100 kms~'Mpc~' and is thought to 
lie in the range 0.5 ^ 1]. Since the first results of the WMAP mission have been made public, we now know 
the constant h with a much smaller uncertainty (Bennett 2003). Cinzano and van der Marel (1994) also showed 
that the total mass of the dark halo is about 1.8 x 10^' h~^MQ. Thus, according to this estimate, the total mass 
of NGC 2974 is 3.0 x lO'^ /2~'Mq, which is obviously smaller than that given by Kim et al., mentioned above. 
The difference appears to be due mainly to differences in the adopted distance and in minor details of the model 
(van der Marel, private communication). Kim (1988) showed that the optically measured heliocentric velocity 
of NGC 2974 is 1924 kms"' according to the optical definition v = cAX/Xq, and that the distance to NGC 2974 
is given in terms of the value of Hubble's constant and is 22.5 h~^Mpc in terms of the Virgocentric flow model 
of Aaronson et al. (1982) [it is 19.25 h'^Mpc based on the Hubble law]. The rotational velocity of NGC 2974 
is about 360 kms~^ (Kim et al. 1988) and its central velocity dispersion is about 222 kms~' (Kim et al. 1988; 
Cinzano and van der Marel 1994). Gauss-Hermite fit parameters of NGC 2974 along the major and minor axes 
have been derived independently by Bender et al. (1994) and Cinzano and van der Marel (1994). 

We use a stellar dynamical galaxy model to calculate Gauss-Hermite fit parameters (a, v, /t3 , 114) as a function 
of the dimensionless position variable x' (measured from the potential centre along the major axis). The physical 
parameters in the models can then be related to the physical dimensions of the observed galaxy with the help of 
the length and velocity scales determined by minimizing a certain function. This function can be constructed in 
the physical domain by use of dimensionless Gauss-Hermite fit parameters {a,v,hi,, /14) , together with the corre- 
sponding observations of an elliptical galaxy, by using a weighted least squares method. An isotropic model can 
then be selected from this class of axisymmetric elliptical galaxy models by analysing the image and ellipticity 
of an elliptical galaxy (in our case NGC 2974). It is known from Fig.|4]in Sect.|5]that the Gauss-Hermite fit pa- 
rameters (a,v,/23,/24) of the model can fit the corresponding observations of the elliptical galaxy NGC 2974 in 
terms of the above least squares method. Finally, the total mass of the best fit model of NGC 2974 can be derived 
from the length and velocity scales obtained by comparing the Gauss-Hermite fit parameters {a,v,hi,,h4.) with 
the corresponding observations of NGC 2974. 

The models we give here are ideal, in that effects from the surrounding environment are not considered. 
The more a galaxy is affected by neighbouring objects, the less well the observations can be expected to be 
fitted by such an ideal model. A galaxy needs to be quite close to us for adequate observations to be available 
to determine such effects. It is not obvious that the best fitting model will be a flattened Jaffe model, but here 
we restrict ourselves to the latter. We first construct the oblate Jaffe model and search for a fit. Obviously, the 
existence of a good fit is not guaranteed but it turns out that one can be found. It is possible that other such 
'suitable' galaxies amenable to this analysis exist but, up to now, no suitable galaxy other than NGC 2974 has 
been found. 

This paper is arranged as follows. We first introduce oblate Jaffe models and Gauss-Hermite fit parameters 
in Sect. 121 Then, in Sect. [3] we mainly analyse relations between the oblate Jaffe model and observed quantities, 
in order to obtain the length and velocity scales of the oblate models. Having derived the finite total mass of the 
models from these scales in Sect. |4l we can then, based on the Gauss-Hermite fit parameters and the ellipticity 
of an elliptical galaxy NGC 2974, find an isotropic oblate Jaffe model in Sect. |5] that is a good fit. The derived 
total mass of the isotropic model only reflects the total luminous mass of NGC 2974. This is because we can 
neglect the influence of dark matter in NGC 2974 on the image, ellipticity and Gauss-Hermite fit parameters of 
the observed galaxy within the radius of 0.5 Re. Our work is summarized in Sect.|6l 
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2 Oblate Jaffe Models and Gauss-Hermite Fit Parameters 



Jiang (2000) defined a class of oblate axisymmetric elliptical systems with the potential-density pairs 




( 



) 



(1) 



^J + {VF+c^ + dY + rj 



and 




M c^drjT^ + r]x^{Y^ + c^d) + (3t + 2rj)c^rjY {Y + dY 



(2) 



where R and z are two of three cylindrical coordinates {R,z,(j>), M, c, d and rj are positive constants, G is the 
gravitational constant, Y = \/z^ + c^ and t = y^R^ + {Y +dY. Here, these systems are obtained by flattening the 
spherical Jaffe (1983) model and they are oblate, and so they are called oblate Jaffe models. Some more general 
extensions of the spherical Jaffe model have been given by Jiang and other coauthors (Jiang and Moss 2002; 
Jiang, Fang and Moss 2002). It is worth mentioning that the spherical Jaffe model is a special case of a family 
of spherical / models (Kuzmin, Veltmann, Tenjes 1986; Dehnen 1993; Saha 1993; Tremaine et al. 1994) which 
have some similar flattening properties (Jiang, Fang, Liu and Moss 2002; Jiang and Ossipkov 2006; Ossipkov 
and Jiang 2007a,b). In the literature, there are also other spherical models (e.g. Veltmann 1961, 1965, 1979, 
1981; Kuzmin and Veltmann 1967a,b, 1973; Ossipkov 1979a,b; Binney and Tremaine 1987) and generalized 
systems with three-dimensional gravitational potentials (e.g. Kutuzov and Ossipkov 1981). 

It is known that the density ^ at large distances decays radially like r"'*, except on the major axis, and as 
on the major axis, where r is one of three spherical coordinates {r,6,^). The central ratio a of isodensity 
contour /?-axis to z-axis (Jiang 2000) is given by 



The total mass of the model is finite and equal to M. 

Jiang (2000) also proved that the form of the potentials allows their two-integral even distribution functions 
to be calculated from the self-consistent density-potential pairs of the above oblate Jaffe model by using a 
modification of Hunter and Qian's algorithm (Hunter and Qian 1993; Qian et al. 1995). Once the two-integral 
even distribution functions have been obtained, the distribution functions /(e,Lj.) can be found by using the 
maximum entropy principle (Dejonghe 1986), where e and are the relative energy and the z-axis angular 
momentum respectively. For oblate stellar models, other distribution functions can also be found depending 
on whether isotropy or anisotropy (Binney and Tremaine 1987) is assumed. Recently, some new formulae for 
distribution functions for both spherical and axisymmetric galaxies have been presented by Jiang and Ossipkov 
(2007a,b,c); but they do not seem applicable to the oblate Jaffe models. 

In order to understand the Gauss-Hermite fit parameters clearly, we first assume that a galaxy is in the 
frame of reference {x,y,z) with mass density given by Eq. Q, and Ox,y,z are its principal axes. Obviously, 
R = Then, without loss of generality, we take the viewing direction to be in the {y,z) plane and project 

the galaxy along the line of sight direction, which is taken as the z'-axis, with an inclination angle (p to the z-axis, 
onto the perpendicular {x' ,y') plane. Thus a new frame of reference (x',y',z') is formed by the sky plane (jc',}'') 
and the line of sight which is the z'-axis; the frames of reference satisfy the following relations 

X = x', y = y' cos (p — z' sin 9, z = y' sin 9 + z' cos (p. (4) 

When the galaxy is observed at an inclination angle (p, the projected surface density Z{x' ,y') at any point {x' ,y') 
on the sky plane is 

L{x',y')=j p([y]^ + [ycos(p-z'sin9]^,ysin9 + z'cos(p)iiz'. (5) 

Let V|| and ay be the line-of-sight velocity and velocity dispersion respectively. The components of velocity 
directed along (x',/) and the line of sight z' can be obtained as follows: 

V,/ = {xvR~yv^)/R, Vyt = (yv^+xv^)cos9/7? + v^sin9, (6) 
V|| =v,t = — (yv^ + ji:v0)sin9//? + Vjcoscp. (7) 



{c + d){\5c"'+ llc^d +\5cd^ + + \Oc^rj + 9cdr j + 3d^rj) 
c^{5c + d){3c + 3d + 2rj) 



(3) 
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The line-of-sight velocity distribution curve of any distribution function can be given by evaluating a triple 
integral of the distribution function /(e,Lj.) with respect to two velocity variables (say, v^/ and v,,/) and one 
position variable (say, z') when £ and L- are expressed as functions of six variables (say, vy, Vyi, v^/, x', y' and 
z') in the phase space. Without loss of generality, we consider the normalized line-of-sight velocity distribution 
curve 

Zo(v||,y,y)=Lo(v||,/,y)/z(y,y) cs) 

with 

Lo(v||,x',y) = j j lfie,L,)dVy,dv,,dz\ (9) 

where the triple integral is to be performed over all the physical values determined by the potential (see, e.g., 
Evans 1994). 

On the other hand, any line-of-sight velocity distribution curve can be quantified by a Gauss-Hermite series 
with its Gauss-Hermite coefficients h„, which can parameterise the velocity distribution curve to carry most of 
information on the shape of the velocity distribution curve (van der Marel and Franx 1993). /o(v|| ,x',_y') is indeed 
approximated by the following Gauss-Hermite series: 

Z(v||,y,y) = 7a-^£/2„M„(w), (10) 



with 



h„ = h„ix',y') = ^^^ I t/v||Zo(v||,x',y)M„H(n = 0,l,2,3,---), (H) 



and 



Uniw) = {2"+'7tnl)-2H„{w)e-- . (12) 

Here w = (v| | — v) /a and H„ (w) (n = 0, 1 , 2, • • • ) are the standard Hermite polynomials. Thus the Gauss-Hermite 
coefficients h„ can be used to explain the properties of the line-of-sight velocity distribution curve of the distribu- 
tion function of stars in any elliptical galaxy, ho represents a Gaussian profile, and odd and even Gauss-Hermite 
coefficients h„ describe asymmetric and symmetric deviations from a Gaussian profile respectively, v is called 
the mean radial or streaming velocity, and a the velocity dispersion (van der Marel and Franx 1993). In fact, 
the quantities v and a are two free parameters and differ from the line-of-sight velocity vy and dispersion C7||. y 
is also a free parameter, in addition to v and a. The parameter y is called the line strength (van der Marel and 
Franx 1993). 

The Gaussian best fit to the velocity distribution curve is such that ho = I and hi = h2 = with an addi- 
tional condition that the fit parameter y is usually around 1 . The Gaussian best fit parameters (v, a, y, h„) can be 
calculated enabling comparison with the corresponding observations of elliptical galaxies. When the Gaussian 
best fit parameters (v, a, fjin) are obtained, improved estimates for the velocity moments can also be given and 
so the Gauss-Hermite coefficients h„ are also termed the Gauss-Hermite moments. Since ho — 1 =h\ = /t2 = 
and yc^ 1 for the above Gaussian best fit, the other four Gaussian best fit parameters (v, a,/23,/24) determine 
principally the velocity distribution curve, and thus they are usually called the Gauss-Hermite fit parameters. In 
fact, once it is assumed that ho = \, we can regard y as a function of the two parameters; v and a, that is, we 
know from ( fTTT) that 

y= V2l_Jv\\loiv\\,x',y')e-'T, (13) 
which gives that < y < a/2 for each point {x',y'). 
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3 Relations between the Model and the Observed Variables 



In order to apply the oblate models to analyse a galaxy, we first have to know the relations between the model 
variables {c,d,rj) of Eq. ([Hi (or (|2]l) and the physical dimensions of the observed galaxy. The stellar dynam- 

1 /2 

ical model takes rj and u = [GM/rj) ' as units of length and velocity, respectively. Thus, in order to relate 
{c,d,rj) in expression ([T| to the physical units of the galaxy, we first have to change {c,d,rj) into dimension- 
less parameters {c,d), assuming that the unit of length is rj [i.e. using a transformation c = crj,d = drj]. Then 
a stellar model with dimensionless parameters {c,d) can be obtained [that is, dimensionless parameters {c,d) 
can be determined from the relation ©I in terms of the image and ellipticity of the galaxy. Finally, the length 
and velocity scales rj and u can be determined by comparing the Gauss-Hermite parameters of the model with 
those of NGC 2974, using a weighted least squares method. This comparison is made by minimizing a function 
constructed below in terms of the weighted least squares method. Thus we can relate the physical parameters 
{c,d, rj) of the models to the physical scales of the galaxy. 

Assume that (a,v,/j3,/j4) are Gauss-Hermite fit parameters, as functions of the dimensionless position vari- 
able x' measured from the potential centre along the major axis for the dimensionless model with variables (c, d) . 
Then {a,v,hj,h4) are also dimensionless. Using the linear relation 

r = rjx' ,v = uv,a = ua. (14) 

we can relate the physical parameters of the models with the physical dimensions of the galaxy. 

The observations corresponding to the Gauss-Hermite fit parameters (a, v, hi,, h^) are functions of the dis- 
tance r from the galactic centre. Generally, 1 arcsec is chosen as the unit of the distance r and 1 kms~' as the 
unit of velocity. The dimensionless model variables a and v are functions of dimensionless position x'. For the 
oblate models the dimensionless coordinate x' is defined as the distance from the potential centre. In this paper, 
it is also assumed to be a position variable from the potential centre along the major axis since we are fitting 
along the major axis only. In order to construct a fit for a galaxy, we have to use the linear transformations (114b 
to convert the dimensionless model variables a, v and x' into dimensional model variables a, v and r, which can 
be compared with the corresponding observations. 

From Eq. (1141) . the Gauss-Hermite coefficients and in the models can be rewritten as 

/J3 = h^{x') = h^{r/rj)M = hA{x') = h^{r/rj)- (15) 
and the dimensional model variables a and v can be regarded as functions of the distance r, i.e., 

V = uv[r/rj) ,a = ua [r/rj) . (16) 

Once the length and velocity scales rj and u have been determined, the physical parameters of the models 
can be found from Eq. (1141) . in order to describe the physical scales of the galaxy. In order to determine the length 
and velocity scales rj and m, we have to use the weighted least squares method to construct a function. Then the 
two scales rj and u can be obtained by minimizing the function in the physical domain. The detailed process of 
obtaining the length and velocity scales rj and u is presented below. 

Observations corresponding to the Gauss-Hermite fit parameters for any galaxy are discrete. Assume that 
there are a sequence of observations [say, a,-, v,-, /la;, h^i] with errors (Sa, , 5v/, 5/i3, , dh/^i), corresponding to the 
Gauss-Hermite fit parameters, at distances r, [i — 1,2, • • • Let /(o) and g{rj,u) be given by the following 
two functions: 

f{rj) = £ (w3,-|/j3(r,yo) -%|' + H'4,i/j4(n/o) -/J4,i') (17) 
(=1 

and 

n 

gin,"-) = ^ {wvi\uv{ri/rj) -Vi\'^ + Wai\ua{ri/rj) ~ a,-p) , (18) 
)=i 

where the weights W3, , W4, , w,,, and Wai are related to the observational errors. Thus the length and velocity 
scales rj and u can be determined by minimizing f{rj) +g{rj,u). 

The scale lengths rj for the arguments v and a are the same as those for the arguments h„. The scale lengths 
are different from the effective radius R^, or core radius R^. Here, the effective radius R^ is defined to be the radius 
of the isophote containing half of the total luminosity and the core radius R^ is where the surface brightness has 
fallen to half of its central value. It is found that the scale length rj is about half that of the core radius [i.e. 
rj^Rc/2] for NGC 2974. 
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Generally, the weights ws, , W4, , Wyi and Wai are chosen as follows. 

W3i = e-45/>3,/A3» ^ = g-45/>4,/''4,„ ^ (19) 

e — 1 e — 1 

WW = , g-^^"', w,,- = , e-'^"-. (20) 

e — 1 e — 1 

We take /j3,„ = /j4m = 0.2 and v„, = Om = 300 kms~\ since the range of the measured velocity dispersions 
is usually 100 - 300 kms"^ (Gerhard 1993a; Bender et al. 1994) and the third-order and fourth-order Gauss- 
Hermite coefficients range from —0.2 to 0.2 (Gerhard 1993a; Surma and Bender 1995) for elliptical galaxies. 
Our choice of the weights (119b and ( 120b is based on the following three reasons. First, the weights should be equal 
to zero when the observational errors attain their corresponding maximal values (0.2 for the Gauss-Hermite 
moments and 300 kms^' for rotation velocities and velocity dispersions), and unity when the observational 
errors are zero. Secondly, the weights should decrease with the observational errors. Thirdly, the weights should 
be continuous functions, for convenience of calculation. 



4 Determination of tlie Total Mass of tlie Model 

It has been shown in Sect. [3] that the model variables can be related to the corresponding observations of the 

1 /2 

galaxy if the length and velocity scales rj and u = [GM/rj) ' are known. In this section we show how to use 
the two scales rj and u to calculate the finite total mass M of the oblate model defined by Eqs. (HJ and 

As mentioned in Sect. [3] the distance r from the galactic centre is generally given in arcsec (Cinzano and 
van der Marel 1994; Bender et al. 1994) and the model variable jc' is dimensionless. For convenience when 
comparing the model variables with the corresponding observations, the two scales rj and u in ( 114b are given 
in arcsec and kms~ ' respectively. Suppose that is the distance from this galaxy to the observer in units of 
Mpc. Assume that p (in arcsec) is the length scale rj of a model. If the Gauss-Hermite fit parameters for the 
model with the two scales rj and u can fit the corresponding observations of a galaxy following minimization of 
a function with ( 117b and dlSb . then p and u can be obtained by comparing the model variables (x', v, a,hj,,h4) 
with the corresponding observed quantities {r,v, a,/i3,/j4) by use of the linear transformation (114b . On the other 
hand, the length and velocity units of the model are rj and (GM/rj)^^^ . Since the two scales rj and u, obtained 
by comparing the Gauss-Hermite fit parameters of the models with the corresponding observations, are the same 
as the length and velocity units of the model respectively, we have 

^^=^"( 180x^3600 ) ^^'^ 

and 

{GM/rjY^^ = u. (22) 

Indeed, (121b is also a transformation from a distance unit of 1 arcsec to a unit of 1 Mpc. Since it has been shown 
by Jiang (2000) that the total mass of the model is finite and equal to M, it is known from Eqs. ( 121b and ( 1221) that 
the total mass M of the model is 

M^r/ u^dn f np \ 



5 Fitting of Oblate Jaffe Models 

In this section we obtain a good fit of an oblate Jaffe model by mimicking the image and ellipticity of the galaxy 
NGC 2974, and comparing the Gauss-Hermite fit parameters of the model with the corresponding observations 
of the galaxy. But, it is finally found that the derived finite total mass of the model is quite close to the total 
luminous mass of this galaxy. 

It has been mentioned in Sect. [3] that the oblate model defined by Eqs. ([T]| and (|2]i with physical parameters 
(c, d,rj) first has to be nondimensionalized to the dimensionless variables c and li, which are given by analysing 
the ellipticity and contours of the image of the galaxy, so that the model can mimic the image of the galaxy. 
NGC 2974 is classified as E4, and the flattening of the contours in B-R images of NGC 2974 (Kim 1989) 
suggests that the inclination angle of the minor axis to the line of sight of NGC 2974 is about 60° (Cinzano 



Fig. 1 The image of NGC 2974 obtained from NASA/IPAC 
Extragalactic Database (NED). 




and van der Marel 1994). Thus we examine oblate models of ellipticity E4 with an inclination angle of 60° to 
describe NGC 2974. On the other hand, it is easy to see from Eq. OJ that the dimensionless model is of type E4 
when c — 8.43 and d=\5 [that is, c = 8.43rj and d = 15rj]. Of course, there are an infinity of solutions {c,d) of 
Eq. (O that make the dimensionless model be of type E4, and {c,d) — (8.43, 15) is just one of these, and is such 
that not only is the ellipticity of the model close enough to that of the galaxy NGC 2974 but also the contours of 
the projected surface density of the model are quite similar to those of the image of the galaxy. That is, it can be 
applied to simulate the image of the galaxy. We find the image of the galaxy from the NASA/IPAC Extragalactic 
Database and plot the contours of the mass density of this model and its edge-on projected surface density. Here, 
Fig. [T] shows the image of NGC 2974, and Fig. |2] and Fig. [3] display the contours of the mass density and its 
edge-on projected surface density, respectively. In Fig. |2] and Fig. [51 successive contour levels differ by factors 
of 0.5, from the potential centre to the outskirts, and the isophotes get more flattened outwards along the major 
axis. This illustates that the isophotes of the model become more disky outwards, as the mass density drops off 
slowly like along the major axis. 

Once the dimensionless model is determined by using the image and ellipticity of the galaxy, the Gauss- 
Hermite fit parameters at any inclination angle for the model can be calculated, by using the maximum entropy 
principle or assuming the isotropy (or anisotropy) of the model (Jiang 2000). Then we use a weighted least 
squares algorithm, as mentioned in Sect. [3] to compare the Gauss-Hermite fit parameters of the model with the 
corresponding observations of NGC 2974. Since the model is oblate, we can assume arbitrarily that the model 
is isotropic. Minimization of the function f{rj) +g(rj,u) shows that an oblate model given by Eqs. ([l]) and (|2]l 
with c = 8.43o and d = 15rj can be used to determin the Gauss-Hermite fit parameters of NGC 2974, with an 
inclination angle of 60°. Further, the length and velocity scales are rj = 0.55 arcsec (or say, p = 0.55 arcsec) 
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Fig. 4 Gauss-Hermite fit parameters for NGC 2974: the error bars represent the observations of Cinzano and van der Marel (1994), 
the broken curves come from the oblate model defined by Eqs. (TJ and (|2} with c = 8.43o and d = 15o, assuming isotropy with 
the choice of rj = 0.55 arcsec and u = 3300 kms^' . The inclination angle of the minor axis to the line of sight is 60°. 



and u = 3300 kms^' respectively. Fig. |4] shows this fit. Our diagnostic or criterion for deciding the quality of 
the fit is that the weighted least squares fit is regarded as good if the mean absolute error values of the velocity 
V, the dispersion a, and the Gauss-Hermite coefficients hj, and /i4, between the model and the observations, 
are less than 15 kms^', 15 kms^', 0.01 and 0.01, respectively. For this model, the mean absolute error values 
between these fit parameters (v, o,h^,h/s) and their corresponding observations are, respectively, 12.75 kms~\ 
12.00 kms"' , 0.005 and 0.004. Therefore this fitting can be regared as good. 

The error bars in Fig.|4]represent the major axis Gauss-Hermite fit parameters of NGC 2974 [see Cinzano and 
van der Marel 1994, figure 1 or figure 5(a)]. These observations are obtained by analysing the galaxy spectrum 
data with the modified Fourier fitting method developed independently by van der Marel and Franx (1993) and 
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Gerhard (1993b). This data analysis is used to determine the derivations of the velocity distribution curve from 
a Gaussian. 

The galaxy spectrum is assumed to be the convolution of a template spectrum and the velocity distribution 
curve. The velocity distribution curve is expanded as the Gauss-Hermite series of the five Gauss-Hermite fit 
parameters, which are defined by dlQI) with (1111) and (1121) . The parameters are determined by 'j^^ fitting in Fourier 
space of the broadened template spectrum to the galaxy spectrum. The finite spectral resolution of the instrument 
needs a high- and a low- wavenumber cut-off. To obtain parameters (7, v, a) from the galaxy spectra, a Gaussian 
velocity profile can be fitted to the spectra. This requires less rebinning of the data in the spatial direction. To get 
the parameters {h^ , h/^), the parametrization (112b can be fitted to the data. 

The galaxy spectra are in fact the data that are produced by the standard reduction of stellar spectra with 
one spectrum of high signal-to-noise ratio for each star. This data reduction is done as follows. First, wavelength 
calibration is done using He-Ar arc lamp spectra, and so is sky subtraction using the available data at the ends 
of the slit. After that, frames should be flat-fielded and cleaned from cosmic rays and CCD defects. Finally, 
the galaxy spectra are rebinned in the spatial direction, with emission lines of 4959 A, 5007 A and 5200 A 
interpolated over, in order to increase the signal-to-noise ratio. 

The stellar spectra are obtained by using a spectrograph equipped with CCD. To do this, one has to consider 
the exposures of the spectrograph, the scale of the rebinned stellar spectra, the sky seeing, the dispersion of the 
grating, the range of the spectra, the width of the slits, and so on. The spectroscopic data of slits along the major 
axis (PA=225°) of NGC 2974 are first taken through two separate 45-min exposures that are reduced separately 
and added later. The stellar spectra are then rebinned at the telescope over 3 pixels in the spatial direction and 
so have a scale of 1.65 arcsec/pixel. The seeing during the observations is about 1.9 arcsec FWHM. The grating 
used has a dispersion of 0.89 A/pixel. The chosen spectral range is from 4750 to 5600 A, centred at about 5175 
A. The slit width is 1.51 arcsec so that the instrument velocity resolution is 45 km/s. 

The broken curves in Fig.|4]are those for the parameters of the oblate model mentioned above. The curve of v 
fits well the velocity of NGC 2974 at radii between 7 to 14 arcsec, although inside a radius of about 7 arcsec, the 
model velocity is lower than that observed, and is higher at distances larger than about 14 arcsec. The curve of a 
lies within the corresponding error bars in the region from 11 to 16 arcsec, but it is lower inside a radius of about 
2 arcsec and higher in the range from 2 to 1 1 arcsec. As the radius increases, the dispersion of the model always 
decreases; but the dispersion of NGC 2974 tends to increase with radius from about 12 arcsec. The graphs of 
the parameters and h/[ are the best fit for the four broken curves in Fig.|4l they are within, or outside but very 
close to their corresponding error bars. In a word, the fit of the four parameters in Fig.|4]is credible, and so it is 
a good example to demonstrate our approach of modelling elliptical galaxies dynamically. 

The Virgocentric flow model of Aaronson et al. (1982) gives do = 22.5 h"' Mpc. Eq. (EB, when p = 0.55 
arcsec and do = 22.5 Mpc, gives rj = 0.06 kpc. Then, by using (123b , the total mass M of the model, given by 
c — 8.43o and d = \5rj, can be calculated as M = 1.52 x lO" /j^'Mq. This mass is very close to the total 
luminous mass of NGC 2974 mentioned in Sect. [T] [i.e., that predicted by Cinzano and van der Marel (1994)]. 
This is consistent with the total mass M of the models only reflecting the total luminous mass. This is partly 
because we choose the model only in terms of the image and ellipticity of the luminosity of the galaxy, and 
partly because the influence of the dark halo in NGC 2974 on the Gauss-Hermite fit parameters of the galaxy 
is negligible within a radius of 0.5 /?e, and the fitting of Gauss-Hermite fit parameters given in Fig. |4]is made 
just within this region where the dynamical influence of the dark halo is negligible (cf. also Cinzano and van der 
Marel 1994). But the radius of the isophote containing half of the total mass of this model is 18.59 arcsec, which 
is less than the effective radius of NGC 2974. This also implies that the dark matter halo extends to the central 
area of the ellipsoidal galaxy NGC 2974. It is worth mentioning that Weijimans et al. (2008) recently found that 
within 5 /?e at least 55 per cent of the total mass is dark for the early-type galaxy NGC 2974. 

Now we consider another model of type E4. When (c, d) = (1 .26, 1 ) [that is, c = 1 .26^ and d = rj'\, it can be 
found from Eq. (|3]l that the dimensionless model is of type E4. Fig.|5]and Fig.|6]show the contours of the density 
of this model, and its edge-on projected surface density, respectively. They are similar to the corresponding 
figures for c — 8.43r7 and d = I5rj (see Fig. |2]and Fig. |3]l, respectively. Similarly, in order to fit the galaxy 
Gauss-Hermite parameters, the length and velocity scales (o,m) can be calculated as follows: rj = 3.2 arcsec 
(or say, p = 3.2 arcsec) and m = 1 1 10 kms^' . Fig.|7]shows a Gauss-Hermite parameter fit for NGC 2974. It can 
thus be found that the total mass M of this model is equal to 1.0 x lO'' /z^'M©. This mass is less than that of 
the model with {c,d) = (8.43, 15), and it is farther away from the total luminous mass of NGC 2974 than that 
of the latter. On the other hand, in Fig. |4] and Fig. |7] the corresponding broken curves of the two models are 
almost similar and their difference is very small. Our criterion of determining a better model is thus as follows: 
the closer the total mass of the model is to the mass of NGC 2974, the better the model is. Therefore, we believe 
that the model with {c,d) = (8.43, 15) is a better fit for NGC 2974 than that with {c,d) = (1.26,1), although the 
velocity v fitting in Fig.|4]is worse than that in Fig.|7] For the model with {c,d) = (1.26, 1), the mean absolute 
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Fig. 5 As Fig.|2]but for c = 1.26o and d = rj. Fig. 6 As Fig.|3]but for c = l.lbrj and d = rj. 



error values between the parameters (v, a^h^^h^) and their corresponding observations can be calculated as 8.75 
kms~\ 14.00 kms~', 0.0085 and 0.004, respectively. Of course, the fitting in Fig.|7]is also a good fit. 

In Fig. |4] and Fig. |7] the systemic velocity of the two models is monotonic while the corresponding galaxy 
observation values flatten out. This is a shortcoming of our models. It might result from our assumption of 
isotropy, which would in turn imply that the galaxy is very anisotropic. 

It is further found from our computations that models of type E4 have similar profiles of the Gauss- Hermite fit 
parameters (v, CT, , ^14) . Generally, the smaller the ratio c /c?, the better are the models of the type corresponding 
to that of the galaxy. Hence, there might be other better fit models for NGC 2974, although we cannot discuss 
these possibilities here. 



6 Conclusions 

Dynamical models are generated from both the gravitational potentials and the distributional functions of the 
stellar systems, based on dynamical principles. Therefore the analysis of the dynamical model described above 
can be summarised into two steps. The first is to search for a dimensionless flattened Jaffe model whose projected 
surface density generated from its gravitational potential can be used to fit the actual surface brightness of the 
galaxy. To do this, the ellipticity of the chosen model must be close enough to that of the image of the galaxy that 
the contours of the model are similar to those of the image of the surface brightness of the galaxy. Of course, the 
contours of the projected surface density of the model must become more disky outwards as its corresponding 
density drops off slowly along the major axis. The second step is to determine the length and velocity scales of 
the dynamical model by comparing the model variables with the corresponding observational scales (i.e. by com- 
paring their Gauss-Hermite fit parameters which are needed for quantifying their velocity distribution curves). 
The Gauss-Hermite fit parameters of the model are obtained from its distributional functions. 

The ellipticity of the model can be determined by the ratio of /?-axis to z-axis extent of the contours of 
the density p{R^,z) near the origin. In order to compare the fit parameters (v,a,/j3,/j4) of the model with the 
corresponding observations, the ellipticity of the model should be the same as that of the B-R images of the 
considered galaxy. The ellipticity of the galaxy being considered can be obtained from observations. By using 
the ellipticity obtained, a model can be chosen to fit the image of the surface brightness of the galaxy being 
considered. For example, for the oblate model, c can be obtained from Eq. OJ if it is assumed that c = crj and 
that d = rj. The model has first to be transformed into a dimensionless model in terms of c, and then the fit 
parameters (v, a,hT,,hn) of the dimensionless model have to be calculated. It is necessary to make an assumption 
about the anisotropy of the models; we choose an isotropic model. 

Comparison with observations requires relations between the observational scales and the model variables. 
The relations are given by the linear transformation (114b . This transformation can be determined if ry and 
(GM/rj)'/^ are known. Therefore, comparing the model variables with observations by use of the linear trans- 
formation, a good pair of values of rj and {GM/rjY^^ can be chosen, such that the model variables fit the 
corresponding observations. 

The total mass of the flattened Jaffe models can be shown to be M. Thus the values of rj and {GM /rj)^^^ 
finally lead to the total mass M of the model. Since the model only reflects the B-R images of the galaxy being 
considered, namely NGC 2974, the total mass M of the model is obviously less than the total mass of this galaxy 
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Fig. 7 As Fig.O but the broken curves are obtained for c = 1.26ry and d = rj, assuming isotropy, with the length and velocity 
scales rj = 3.2 arcsec and = 1110 kms^' . 



and so it reflects the total luminous mass of this galaxy. Additionally, the dynamical influence of dark matter in 
NGC 2974 on the Gauss-Hermite fit parameters of the galaxy observed is negligible inside a radius of 0.5 R^. 

For a very flattened model, the two values of rj and {GM/rjY^^ can be also determined by comparing the 
circular velocities generated from the potential of the model with data of the rotation velocities of the very 
flattened galaxies observed. This is the customary approach of fitting the projected surface density of the model 
to the observed surface brightness of the galaxies. This result must be the same as that obtained by comparing 
the Gauss-Hermite fit parameters, if the model is a very good fit. 
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The resulting fit in Fig.|4]is satisfactory, and so it is a good example to illustrate our metliod for constructing 
such models. And, up to now, no other plausible alternatives have been found for choosing a suitable flattened 
Jaffe model for the elliptical galaxy NGC 2974. 
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